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ABSTRACT 

We estimate cluster ages from lithium depletion in five pre-main-sequence groups found within lOOpc of the 
Sun: TW Hydrae Association, r\ Chamaeleontis Cluster, (3 Pictoris Moving Group, Tucanae-Horologium As- 
sociation and AB Doradus Moving Group. We determine surface gravities, effective temperatures and lithium 
abundances for over 900 spectra through least squares fitting to model-atmosphere spectra. For each group, we 
compare the dependence of lithium abundance on temperature with isochrones from pre-main-sequence evo- 
lutionary tracks to obtain model dependent ages. We find that the 77 Chamaelontis Cluster and the TW Hydrae 
Association are the youngest, with ages of 12±6Myr and 12±8Myr, respectively, followed by the j3 Pictoris 
Moving Group at 21±9Myr, the Tucanae-Horologium Association at 27±1 1 Myr, and the AB Doradus Mov- 
ing Group at an age of at least 45 Myr (where we can only set a lower limit since the models - unlike real stars - 
do not show much lithium depletion beyond this age). Here, the ordering is robust, but the precise ages depend 
on our choice of both atmospheric and evolutionary models. As a result, while our ages are consistent with 
estimates based on Hertzsprung-Russell isochrone fitting and dynamical expansion, they are not yet more pre- 
cise. Our observations do show that with improved models, much stronger constraints should be feasible: the 
intrinsic uncertainties, as measured from the scatter between measurements from different spectra of the same 
star, are very low: around 10 K in effective temperature, 0.05 dex in surface gravity, and 0.03 dex in lithium 
abundance. 

Subject headings: open clusters and associations: individual (rj Chamaelontis, TW Hydrae, (3 Pictoris, Tucanae- 
Horologium, AB Doradus) — stars: abundances — line: profiles — stars: pre-main- 
sequence 



1. INTRODUCTION 

Triggered largely by the discovery of young stars in the 
ROSAT X-ray Satellite All-Sky Survey, over the last decade 
several nearby pre-main-sequenc e (PMS) star groups hav e 
been identified (for a review, see IZuckerman & S ong 2004). 
Ranging in age from roughly 6 Myr to ~ 100 Myr, five sep- 
arate groups can be distinguished: TW Hydrae Association 
(TWA), -q Chamaeleontis cluster (77 Cha), (3 Pictoris Mov- 
ing Group (BPMG), Tucanae-Horologium association (TU- 
CHOR) and AB Doradus Moving Group (ABD). The com- 
mon space motions and localized sky positions suggest that 
these groups are likely connected to the Sco-Cen star form- 
ing region located ~1 00p c away in the sout hern hemisphere. 
iMamajek et~ai1(ll999l) and lSong et ai1(l2003l) have traced back 
the space motion of members in BPMG, TWA and 77 Cha and 
argue that the groups are related to a star formation burst in 
the Sco-Cen region as a result of the passing of the Carina arm 
~60Myr ago. 

These groups, because of their close vicinity, are excel- 
lent laboratories for studying star and planet formation. Well 
constrained ages are necessary to make conclusions about 
timescales of, e.g., disk dissipation and planet formation. 
Already, observations from the same sample presented in 
this paper have revealed that ac cretion disks can last up to 
~10Myr, but beyond this is rare (ETayawardh ana et aTl F2006). 

Previously, ages have been derived from Hertzsprung- 
Russell (HR) diagram fitting, group dynamics and lithium 
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abundance measurements. iLuhman & Steeg hs (2004]) pro- 
vide an HR diagram isochrone age for 77 Cha of 6^ Myr de- 
rived fro m the evo lutio nary models of iBaraffe et all (119981: 
hereafter lBCAH98l) and lPalla & Stahleri d 19991) . which agrees 
we ll with the dynamical expansion age of 6.7 Myr determined 
by iJilinski et al.l (120051) . The dynamical age of TWA has 
been harder to determine because of inconsistent space mo- 
tions among its more than 30 members. An inferred age of 
8 . 3 ± . 8 Myr is g iven to TWA based on t he dynamic al motion 
of four members (Ide la Reza et al.ll2006l) . However, the likely 
complex dynamical evolution of TWA has led to sev eral plau- 
sible evolutionary scenarios. Mak arov et al.l (120051) attribute 
this complex evolu t ion to a chance encounter with Vega, while 
lLawson & Crausd (|2005) suggest TWA is composed of two 
separate groups, based on bimodal rotation period distribu- 
tions with distin ctly separate ages of ~ 10 M yr and ~ 17 Myr. 
More recently, Barrado y Navascues (2006) finds a conserva- 
tive age of lO^j Myr by comparing ages from HR diagram 
isochrone comparisons (from BCAH98) and lithium abun- 
dances. 

The slightly older group BPMG has an estimated age 
of 12 ^ Myr bas ed on HR diagram isochr one comparisons 
(from BCAH98) and lithium abundances (Zuckerm an et al.1 
1200 lb . with three dimensional motions that a re consistent 
with a dynamical expansion age of 11.5 Myr dOrtega et al.l 
2002). Feigels on et al.1 ((2006) independently derived an age 
of 13^ Myr for the recently confirmed wide binary system of 
5 1 Eri and GJ 3305, part of BPMG. 

Known to be older than BPMG, but younger than the 
Pleiades, TUCHOR has an age of 20-40 Myr based on Ha 
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measurements, X-ray luminosity, rotation and lithium abun- 
dances in co mparison to other young clusters like TWA and 
the Pl eiades dZuckerman & Webbll2000HStelzer & Neuhauseri 
l2000h . 

Perhaps, the most de bated age is that of the ABD group. 
IZuckerman et al.l (120041) derived an age of 50 Myr by compar- 
ing the Ha emission strength of ABD members to members 
of the younger TUCHOR association, in addition to fitting its 
three M-type members to HR diagram isochrones. In con- 
trast, |Ohmanetal] ([2005b compared HR isochrones of ABD 
members to those of two well-observed clusters with ages of 
50 and 125 Myr and suggested that the ABD group is coeval 
with the Pleiades at an a ge of 100- 125 Myr. The latter age 
is strongly supported by [Ortega et al.l d2007l) . who compute 
full 3D galactic orbits of ABD and the Pleiades cluster and 
show the dynamics of the two groups can be traced back to a 
common origin of 1 19 ± 20 Myr ago. 

A relatively new approach to age estimates is to use 
the evolution of the lithium abundanc e for low-mass, par- 
tially and fully convectiv e PMS stars (Bil dsten et al.l 119971 
iJeffries & Oliveiral [2005). The initiation and duration of 
lithium depletion in PMS stars is dependent on mass and is 
very sensitive to the central temperature. Lithium is con- 
verted into helium in p. a reactions in cores of low-mass stars 
when the temperature reaches 2.5 x 10 6 K. The lower the stel- 
lar mass, the longer the time it takes to reach this critical tem- 
perature. For example, a 0.6 M Q star begins to burn lithium 
at an age of 3 Myr, while a lower mass star at 0. 1 M Q begins 
to burn lithium at an age of 40 Myr. Stars with M < 0.06 M 
never reach this typical temperature, while stars with 0.6 M Q 
< M < 1.2M© burn lithium for a short period (l-2Myr) un- 
til a radiative core develops, and more massive pre-main se- 
quence stars do not destroy lithium in their envelope at all. 
The result of these processes is a dip in the lithium abundance 
as a function of luminosity (and consequently, a function of 
effective temperature), only affecting stars with spectral types 
late than F5. As a group ages, this dip becomes deeper and 
widens on the cool end as progressively cooler stars reach the 
critical core temperature. 

The cool end of this dip has been used to date coeval groups 
containing late M-dwarf stars by identifying the lithium de- 
pletion boundary (LDB). The LDB marks the luminosity 
above which all stars will have depleted their lithium. The 
lithium is very quickly depleted in these low-mass stars, so 
the LDB marks a sharp jump from initial to near depleted 
lithium abundances. As the temperature in the cores of PMS 
stars increases in time, the LDB will shift to cooler tempera- 
tures as a cluster ages. LDB ages have been determined for 
the Pleiades (125 ± 8 Myr), a Per (90 ± 10 Myr), IC2391 
(53 ± 5 Myr), and NGC2547 (35 ± 4 Mvf) (IStauffer et ail 
19981: iBarrado v Navascues et aj.| Il999t IStauffer et al.lll999b 
Barrado v Navascues et al.ll2004t Ijeffries & Oliveirall2005l) . 

In this paper, we fit more than 900 multi-epoch, high- 
resolution spectra, introduced in $2] of 121 low-mass PMS 
stars to synthetic spectra created from PHOENIX model at- 
mospheres (see §|3}. We begin by identifying the model with 
the best-fitting surface gravity, logg, and effective tempera- 
ture, r e ff for each observed spectrum, as described in §|4] In 
$5] we use these best-fit T e ff and logg to find the lithium abun- 
dance by fitting model spectra to both the 6 104 A and 6 708 A 
lithium features in the observations, but now using lithium 
abundance as a free parameter. For comparison, we also mea- 
sure the equivalent width (EW) of the 6 708 A lithium doublet 



which allows us to chronologically order the groups based 
solely on empirical measurements. We proceed by comparing 
the measured lithium abundance distribution to th at predicted 
by PM S evolutionary models of IBCAH98I and ISiess et all 
(2000) and estimate model dependent ages for each of the 
PMS groups. 

2. OBSERVATIONS 

High-resolution spectra were obtained during six separate 
observation runs, on a total of 19 nights, between Decem- 
ber 2004 and April 2006, utilizing the MIKE spectrograph at 
the Magellan-Clay 6.5 m telescope on Las Campanas, Chile. 
MIKE is a double echelle instrument, covering two sepa- 
rate wavelength regions. For this study, the red region from 
4 900 A to 9 300 A is used. The raw data were bias subtracted 
and flat-fielded, and before extraction, the scattered back- 
ground in the spectrograph was subtracted by fitting splines 
to interorder pixels. The spatial direction of the projected slit 
produced by MIKE is wavelength dependent and not aligned 
with the CCD. We therefore extracted spectra using a custom 
procedure develped in ESO-MIDAS, that takes into account 
the tilt and optimally extracts the spectrum by iteratively es- 
timating the slit illumination function. For wavelength cali- 
bration, exposures of a Thorium-Argon lamp were used, as 
well as observed telluric absorption lines. A more detailed 
account of the reductions and a log of the observations will 
appear in a forthcoming paper (A. Brandeker et al. 2008, in 
prep.). Multiple spectra for many of the targets in our sam- 
ple were taken in order to search for multiplicity and perform 
variability studies. 

With no binning and using the 0."35 slit, the spectra in this 
study have a resolution of R^60000. The pixel scale was 
0."13pix _1 in the spatial direction, and about 24mApix _1 at 
6 500 A in the spectral direction. Integration times were cho- 
sen so that the signal-to-noise ratio (S/N) >70 per spectral 
resolution element at 6 500 A, except for the brightest stars 
where this would have implied an exposure shorter than 120 s. 
In those cases, we used the longest exposure time shorter than 
120 s that did not saturate the detector, giving (S/N) = 70-500, 
depending on seeing. 

For this study, we only use objects with spectral types later 
than F5, and earlier than M5. Hotter stars do not show 
depleted lithium abundances in their atmospheres, because 
their cores are already radiative when lithium burning starts. 
Cooler stars were not included in the survey because they 
were too faint. We also exclude obvious spectroscopic bi- 
naries. In total we have a sample of 121 stars, with 11 from 
r?Cha, 32 from TWA, 23 objects from BPMG, 35 from TU- 
CHOR, and 22 from ABD. For comparison, we also analysed 
20 radial-velocity standard stars from the field. 

3. MODELS 

We compute synthetic spectra usin g version 14 of the 
PHOE NIX model atmosphere package (Hauschi ldt & Baron] 

19991). The mode l atm ospheres are described in 

Kucinskas etal] d2005l 120061) . In total, we have 240 
model spectra covering a spectral region from 3 200 A to 
10 000 A with 0.03 A spectral resolution, with temperatures 
(r eff ) ranging from 2 500 K to 6 500 K in steps of 100 K, and 
surface gravities (logg) ranging from 3.0 to 6.0 in steps of 
0.5 (i.e. gravities ranging from 10 3 cms" 2 to 10 6 cms~ 2 ). 
For each of these temperature and gravity combinations, we 
have additional models with different lithium abundances 
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(A^Li) ranging from logiVu = 0.0 to 4.0 in steps of 0.5 (where 
the normalization is such that logA^j = 12), which cover the 
6 104 A and 6 708 A lithium absorption lines. All models 
were calculated at solar metallicity. 

The atmospheric models are calculated under the as- 
sumptions of local thermodynamic equilibrium (LTE). 
ICarlsson et al.l (fl994) have systematically analyzed the effects 
of non-LTE on the formation of the Li I line in cool stars. They 
show that that the non-LTE effects can lead to discrepancies in 
lithium abundances measured from the 6 708 A lithium dou- 
blet from about -0.3 to 0.3 dex depending on the depletion and 
temperature. We attempt to minimize the exclusion of non- 
LTE effects by also fitting the 6 104 A lithium feature whose 
non-LTE corrections are opposite to those of the 6 708 A line 
(ICarlsson et al.lll994l) . 

In addition, we do not consider the effects of chromospheric 
and/or magnetic activity. To compensate for the former, we 
exclude from our fitting method emission lines due to chro- 
mospheric activity that are frequently seen in cl assical T Tauri 
and post T Tauri stars (described further in 34. 11 1. We also note 
here that very strong lines — which are completely dominated 
by pressure broadening, a relatively poorly understood pro- 
cess — require special treatment to reproduce the line profile. 
Although the models include attempts to do this, it does not 
always work satisfactory. The resulting mismatches primarily 
affects lines from neutral metals that are highly abundant and 
very optically thick (such as K I and Na I). For less abundant 
neutral metals, such as Li I, the effects should be minimal. 

4. SURFACE GRAVITIES AND EFFECTIVE TEMPERATURES 

We derive surface graviti es and effective t emper atures, fol- 
lowing the methodology of Moha nty et al.l (2004), by fitting 
theoretical spectra to our observed spectra in small, 80 A 
wide, spectral regions that contain absorption lines or molec- 
ular bands that are highly sensitive to surface gravity and/or 
effective temperature. 

4.1. Modifications to the data 

To ensure our measurements would not be influenced by tel- 
luric absorption lines, or by stellar emission related to activ- 
ity, we made some alterations to the observed spectra before 
fitting. For the telluric lines, we first used multi-epoch ob- 
servations of one of our stars (J3 Pic) to identify the stronger 
atmospheric lines and determine their depths relative to the 
continuum. Next, for each spectrum, we marked a fixed spec- 
tral width of data points around the center of each telluric line 
as contaminated; these regions are ignored in all fits done be- 
low. We used widths of the marked regions that depended on 
the strength of the telluric line (as measured for /3Pic): a re- 
gion of 0.2 A was removed for relatively weak lines, which 
had a depth less than 80% of the continuum. Regions of 0.6 A 
were removed for line depths between 80% and 88%, while 
for the strongest, deeply saturated absorption lines, 1 .0 A wide 
regions were removed. 

As our sample consists of young stellar groups, many of the 
stars show signatures of activity, with relatively narrow chro- 
mospheric emission lines found in many of our targets, and 
strong and wide emission lines ind i cative of accretion i n a few 
others dJayawardhana et all l2006t IScholz et all 120071) . The 
chromospheric activity causes emission features at a number 
of wavelengths, with the strengths and locat ions depending 
on spectral types and age (Eisloff el et alJfl990l) . To mitigate 
the effects of activity, we exclude 2.4 A regions around lines 



TABLE 1 

Spectral Regions Fitted 



Line ID 


A (A) 


AA (A) 


T(K) 


Sens. 


Na P 


5 893 


5 850-5 930 


2 500^1 500 


T, logg 


FeP 


5 893, 5 898 








TiO a 


5 847-6 058 






T 


VO b 


7 851-7 973 


7 900-7 980 


4 000-6 500 


T 


CN a,b 


7 916,7 941,7963 








Na P 


8183,8 195 


8150-8 230 


2 500-3 000 


l°gg > T 


TiO ab 


8 432, 8442, 8 452 


8400-8 480 


2 500-6 500 


T 


FeP 


8 440, 8 468 








Ca IP " 


8498, 8 542 


8485-8565 


2 500-6 500 


T 


VO b 


8521, 8 538 









NOTE. — A list of the spectral regions selected for our fitting 
method as described in 34.2I & 34.31 The first and second columns 
identify and locate spectral features within the chosen spectral 
range that show strong sensitivity to temperature and/or surface 
gravity. The third column lists the 80 A region that was selected in 
our fitting procedure. The fourth column indicates what range in 
effective temperature the selected region shows strong sensitivity 
to varying parameters and the fifth column indicates whether the 
r egion is more sensitive to effectiv e temperature or su r face g ravity. 
a lTorres-Dodgen & Weaver! dl993h b lKirkpatrick et al] dl99lh 

that frequently are seen in emission from classical T Tauri and 
post-T Tauri stars dStahler & Pallal 120051) : He I at 5 875.6 A, 
Na I at 5 890.0 and 5 895.9 A, O I at 8446.5 A, and Ca II at 
8 498.0,8 542.1 and 8 662.1 A. 

Some stars also show signs of accretion, which leads 
to strong, wide emission features, potentially leading to 
less accurate fits. Since only three out of eleven ?yCha 
and two out of 32 TWA members show signs of accretion 
dJavawardhana et al.l2006l) . this will affect relatively few stars 
in the sample and in practice we do not find any effect of ac- 
cretion on the uncertainties involved in constraini ng r e ft, logg 
and logA^Li, except for the EW of one star (see 35. jj . There 
are no known accretors in any of the older PMS groups inves- 
tigated in this study. 

4.2. Selection of spectral regions 

We locate spectral regions that ar e strongly sensitive to sur- 
face gr avity and temperature. BothlTorres-Dodgen & Weaver! 
d!993l) and iKirkpatrick et al.l d 199 lb have compiled lists of 
prominent absorption features useful for spectral classifica- 
tion of low-mass stars. For all, we checked whether our model 
spectra (see $3) indeed showed strong temperature and/or 
gravity sensitivity, and selected only those that did. 

In general, the dependence of these features on our param- 
eters is somewhat degenerate, as different T & s and logg com- 
binations ca n fit a given absorption feature equally well. As 
discussed in Mohan ty et al.l ([2004), features that show con- 
trasting dependence on gravity and temperature are needed to 
constrain the parameter space. Thus, we choose combinations 
of molecular absorption features like titanium oxide (TiO) and 
vanadium oxide (VO) with lines from neutral alkali elements, 
like K I, Na I and Mg I. This works well because, for neutral 
alkali elements, an increasing temperature can be compen- 
sated by increasing gravity, while the molecular absorption 
bands show less correlation between gravity and temperature. 
For instance, TiO bands in the vicinity of 7 100 A are much 
more sensitive to temperature than to gravity, while the triple 
headed TiO bands in the range of 8 440 A are more sensitive 
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FIG. 1. — Spectral fits (solid red lines) to r]Cha 10 (black dots), a K7 
star, in the regions listed in Table[T] Quality of fit (x 2 ) contours at 68.3,%, 
95.4 % and 99.99 % confidence for each spectral window are shown in the 
bottom panel with each level given by x^ in X [1 +(2.30,6.17, 18.4)/JV dof ], 
where iVdof is the degrees of freedom in the fit and X^ m i- s me best-fit value. 
The summed \ 2 values for all spectral windows are represented by the tight 
solid black contours. 



to gravity than to temperature. Examples of this dependence 
can be seen in the \ 2 contours shown in the bottom panels 
of Figs. Q] and [2] Of course, these dependencies change with 
temperature, and hence different sets of regions are best for 
different spectral types. In principle, one could choose to fit 
different regions in different temperature ranges, but 



this risks (borne out in practise) that at the borders there 
are false jumps in parameters due to systematic problems 
with the models. Thus, we fit all regions for all stars, which 
gives very tight constraints on T e s and logg for all spectral 
types (and more smoothly varying systematic offsets; see be- 
low). After investigat ing spectral fits of twelve indi vidual 
regions collected from iTorres-Dodgen & Weaver! ( 1 19931) and 
KirkpatricTetd] (1199C we settled on the five separate spec- 
tral regions listed in Table Q] Also listed is the effective tem- 
perature range where the features provide strong constraints, 
as well as whether these constraints are predominantly on T e g 
and/or logg. 

We excluded regions which either tended to overestimate 
\ogg or provided T e ffflogg values extremely offset from the 
majority of spectral regions. For instance, for many of the 
neutral metal lines, particularly K I at 7 665 and 7 699 A, the 
best fits are at surface gravities much higher than are expected 
forlate-type PMS stars. Since these lines are sensitive to grav- 
ity and generally fit the spectrum very well, the entire fit was 
very sensitive to the K I line. Due to this inconsistency be- 
tween the line strengths in the atmospheric models and pre- 
dicted surface gravity, this region and several others were not 
used. We note that we use other neutral metal lines such as 
Na I that could be flawed similarly. In practise though, these 
produced fits that were consistent with the other spectral re- 
gions investigated. 

4.3. Fitting Methods 

All spectra in our study with spectral types later than F5 
(corresponding to about 6 500 K, the highest temperature for 
which we have models) were fitted to the synthetic spectra 
for each spectral window listed in Table Q] Prior to the fit, 
the synthetic spectra were convolved with a gaussian filter 
to match the observed resolution. The fit used a variant of 
the br oadening-function formalism introduced by iRucinskil 
(2002). In this formalism, a least-squares fit is made of the 
observed spectrum to a set of reference spectra that differ only 
in their velocity offset. This way, rotational broadening is au- 
to matically a ccounted for. The difference with the formalism 
of IRucinskil is that the sum of the model spectra at various 
velocities is also multiplied with a polynomial, to account for 
not only the normalization, but also for small errors in the flux 
calibration. For our small wavelength regions, we found that 
a third-degree polynomial sufficed. 

We note that our fitting method effectively introduces a rel- 
atively large number of parameters that are not of direct physi- 
cal interest, viz., the line shape. Since these parameters might 
be covariant with some of the effects of temperature and/or 
gravity changes, the constraints on T e s and logg we find are 
thus not as strong as might be possible if, e.g., instead we 
broadened the synthetic spectra with an analytic broadening 
model (with only the projected rotational velocity as parame- 
ter). In practice, however, the covariance is small, and our er- 
rors are clearly dominated by systematic mismatches between 
the models and the observations. 

With the above least-squares fitting method, we fitted each 
of the five regions from each observed spectrum to 77 of the 
synthetic spectra, covering all seven values of logg (from 3.0 
to 6.0 in steps of 0.5), and eleven temperatures in a range of 
1 000 K around the effective temperature inferred from the ob- 
ject's spectral type (as liste d in Tab les [2H7J we used the tem- 
perature scale for dwarfs of Luhman & Rieke 1998). For new 
companions with unknown spectral types, initial estimates for 
the temperatures were found through trial and error. 
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FIG. 2.— As Fig. 1, but for HIP 32235 (black dots), a G6V-type TUCHOR 
member. 



To determine combined constraints on 7^ and logg for a 
given spectrum, we add up the \ 2 values derived for the five 
regions. The results are similar to fitting the five regions to- 
gether, except that our method leaves greater freedom for vari- 
ations in continuum (which might be expected) or line shape 
(which would not). The difference with fitting the entire spec- 
trum is that we have effectively ignored the parts of the spec- 
trum that are either insensitive to T e g and logg or where the 
models are shown to fit poorly to the observations. 

For our final estimates of the best- fit values of T e s and logg 



for each spectrum, we interpolate in the 7 by 1 1 grid by de- 
termining the minimum of a two-dimensional parabola (of the 
form a + bx+cy+dx 2 + exy+fy 2 ) fit to the sixteen grid points 
with the lowest values of x 2 - 

4.4. Examples 

In Fig. [T]and |2] we show our fits for one spectrum of each 
of the stars ^Cha 10 and TUCHOR member HIP 33235. The 
gaps in the data in the panels centered on 5 890 and 8 525 A 
are due to the removal of emission lines as discussed in 34.11 
Similarly, on close inspection, one sees that multiple telluric 
lines have been removed from the 8 190 A region. 

Overall, the broadened models reproduce the observed 
spectra well. In detail, however, there are clear inconsisten- 
cies between the observed and synthetic spectra, apparent in 
many of the panels. In particular, a number of lines appear 
to be absent in the models, or are clearly too weak. We also 
found examples of the reverse in some other wavelength re- 
gions. 

In the bottom panels of Figs. Q] and [2j the resulting 
68.3%, 95.4% and 99.99% confidence contours in the T eB - 
lo gj? parameter spac e are shown, with levels set according 
to [Press etal.1 (fl992h at xL x [1 + (2.30,6.17, 18.4)/AW], 
where A^of is the degrees of freedom in the fit and Xmin i s 
the best-fit value inferred from the parabolic fit. Contours 
are shown for each spectral window as well as for the total 
X 2 values. For the different spectral regions, one notices the 
different covariances in T e g and logg, and how, by using a 
number of these sensitive regions, it is possible to constrain 
gravity and temperature precisely. Indeed, the contours for the 
summed x 2 distribution (solid black contours) are extremely 
tight. We will see below that different spectra of the same star 
lead to similarly small scatter in the inferred temperature and 
gravity. 

One also notices, in particular in Fig. [2] that the contours 
for the different regions are statistically inconsistent with each 
other, with differences of several 100K in T e s and up to 1 dex 
in logg. These differences likely reflect systematic uncertain- 
ties in the models, similar to what we find for the resulting 
best-fit average values below. 

4.5. Results 

The best-fit model T e g and logg with corresponding errors, 
for each star, including field stars, can be found in Tables 12^7] 
We consider two wa ys of estimating the associated uncertain- 
ties. First, we follow [Press et al.l d 19921) and use the curvature 
of the best-fit two-dimensional parabola to the x 2 values to 
find regions which are enclosed within a level of 68% con- 
fidence. Second, we consider the error in the mean between 
the results from different spectra of the same source taken at 
different epochs (e.g. the standard deviation divided by the 
square root of the number of observations taken). 

For the temperatures, our statistical uncertainties are typ- 
ically around 9K, while the scatter derived from multiple 
observations of the same object is on average 1 1 K. For the 
surface gravities, our statistical uncertainties numbers are 
0.02 dex on average, and the scatter derived from the error in 
the mean from multiple observations is about 0.03 dex. The 
above suggests the true intrinsic uncertainties (dim) in our 
temperature and gravity measurements are very small, about 
10 K and 0.05 dex for a single observation. We will see be- 
low, however, that systematic mismatches as a function of, 
e.g., spectral type, are much larger. 
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FIG . 3 . — Surface gravity, log g, as a function of effective tempera ture, r e (f, 
for all five PMS groups, as well as field stars. PMS tracks from BCAH98 are 
overdrawn, with, from top to bottom, ages of 4, 8, 15, 25, 35, 45 and 150 Myr. 
We represent the external error in the plotted parameters for the ensemble in 
the bottom-left corner. It is evident that there are systematic problems with 
the values of log g derived from the models, particularly for 3 500 < T e f{ < 
4100K. However, the relative ordering is clear: stars that are older, such as 
ABD and TUCHOR members, have higher surface gravities than younger 
stars. The ultrafast rotators and accretors in our sample are marked with 
circles and squares, respectively. 



In Fig. [3] we show the distribution of logg as a func- 
tion of T e ff for each group as well as for the field stars ob- 
served. In addition, we draw isochrones from the IBCAH981 
PMS evolutionary models. One immediately sees there is a 
clear systematic problem in determining logg for stars with 
3500 < T e ff < 4 100K: logg increases with temperature from 
3 000 K to 4 000 K, but at around 4 000 K, it becomes almost 
1 dex smaller, an unrealistic physical trend. The systematic 
differences in logg can result from a few effects. If the reso- 
lution of the model spectra is not much better than the obser- 
vations before smoothing, line depths can be systematically 
off. In addition, stellar activity also introduces systematic er- 
rors, as has been found for young M-dwarf s, where the chro- 
mosp here feeds back into the photosphere (Fuhrmeist er et al.l 
2005). This feedback is not incorporated into the model, and 
may lead to systematic errors of about 0.5 dex in logg. 

To evaluate the dependence of the fitted T e s on logg, we 
re-fit our spectra three more times using the same method 
described in $4.11 except with \ogg fixed to 4.0, 4.5 and 
5.0. The change in T e g across the surface gravity space is 
small in an absolute sense. Going from log g fixed at 4.0 to 
4.5, the average change in temperature within our sample is 
AT e ff = 71 K. It is slightly higher going from 4.5 to 5.0, with 
ATeff = 124K. Changes from logg = 4.0 to 5.0, yield abso- 
lute average changes in T e s of 162 K. 

Thus, the external errors related to the models are much 
higher than the internal errors. Tables 2 to 7 list the best-fit 
surface gravities and temperatures with quoted errors repre- 
senting the error in the mean of multi-epoch observations of 
the object. More conservatively, we consider our external er- 
rors to be 150 K in T e ff and 0.5 dex in logg. 

In Fig. [4] we compare our computed temperatures to those 
obtained by converting spectral types to effective temper- 
atures using the t emper ature scale for dwarf-type stars of 
iLuhman & Riekd (fl998). Generally, the two scales are con- 
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FIG. 4. — C omparison of the temperatures obtained from spectral fit- 
ting (see j]4.3t with those derived from spectral types using the scale of 
Luhman & Rieke ( 1998). We represent the external error in the plotted pa- 
rameters for the ensemble in the top left corner. The ultrafast rotators and 
accretors in our sample are marked with circles and squares, respectively. 

sistent within the external errors just discussed, with the 
largest deviations show by the accretors (marked by squares) 
and ultrafast rotators (marked by circles, see also jj5.41 i. In ad- 
dition, most of the objects with r e ff > 5000 K, while individu- 
ally consistent within the uncertainties, appear systematically 
to have fitted temperatures slightly greater than those inferred 
from spectral type, suggesting a small systematic error. 

5. LITHIUM ABUNDANCES 

We analyse lithium abundances using two independent 
methods. In one meth od, which follows previous work 
dJeffries & Oliveirall2005l) . we measure the equivalent width 
(EW) of the 6 708 A lithium absorption doublet for each spec- 
trum. We use this measurement and the spectral types quoted 
in the literature to chronologically order the groups. The other 
method uses the best-fit model calculated in $4] to fit model 
spectra for various lithium abundances to regions surround- 
ing the 6 104 and 6 708 A lithium absorption features. As with 
our EW measurements, we chronologically order the groups. 
Further, w e c ompare the obse rvations with PMS models of 
IBCAH98I and lSiess etafldlOOOh to see if the isochrone ages 
based on lithium abundances are consistent with other age de- 
terminations of the five nearby PMS groups in this study 

5.1. Lithium equivalent widths 

First, we use the EW of the 6 708 A lithium doublet feature 
in each observed spectrum as an empirical measurement of the 
lithium abundance. To measure it, we interactively chose the 
edges of the lithium feature and the boundaries of sufficiently 
large regions around it to define the stellar continuum. We 
chose to consistently reject the lowest 15% of the flux points 
in the continuum region (corresponding to the deepest absorp- 
tion features), so that we do not underestimate the continuum 
flux. We do not apply any correction for the 5 identified ac- 
cretors in our sample, but find that veiling due to accretion 
can reduce the measured EW, as discussed briefly below. 

The resulting mean EWs, averaged over all available spec- 
tra, are listed in Tables [2] to [6] with uncertainties being the 
error in the mean between multiple spectra. Fig. [5] displays 
the resulting EWs as a function of spectral type. From this 



7 



is 
w 



0.6 



0.4 



0.2 



i i i i i i i i i I i i i i i i i i i I i i i i i i i i i I i i i i i i i i i 

■ O TWA 
_ i i) Cha 




I I I 



F0 GO KO MO 

Spectral Type 

FIG. 5. — Equivalent width (EW) of the 6708 A lithium line as a function 
of spectral type. Without depletion, the EW would increase with decreasing 
temperature, as is seen with the youngest two groups, TWA and r\ Cha. But in 
older groups lithium is being depleted, and the order, from youngest to oldest, 
is evident: BPMG, TUCHOR, and ABD. We identify the ultrafast rotators 
and accretors in our sample with black circles and squares, respectively. 

purely empirical figure, the chronological order of the groups 
is evident: from oldest to youngest, they are ABD, TUCHOR, 
BPMG, TWA and r \ Cha. The same ordering was found by 
Zuckerman & Song ( 2004) from lithium EWs for a smaller 
sample of stars. Although some TWA members appear to be 
as young as those in 77 Cha, it is quite clear that no mem- 
bers are older than BPM G, contrary to the suggestion by 
lLawson & Crausd (|2005'). (Note that we implicitly assume 
here the initial lithium abundance was the same for all groups. 
We return to this below.) 

We identify the ultrafast rotators and accretors in our sam- 
ple with black circles and squares, respectively. The EWs 
of the rotators in BPMG (red squares) and TUCHOR (green 
stars) stand out above the general trend for each group; we 
will return to this in 35.41 We also note that one accretor in 
77 Cha, (77 Cha 13) has a lithium EW that is lower than most 
objects at the same temperature (~ 3100K) and age, but we 
will find below a lithium abundance near initial from the fits 
to the spectra of this star (see 35.21 ). This suggests the low 
estimate of the EW is related to the fact that it is accret- 
ing, either by its effect on the continuum level or alterna- 
tively, by us underestimating the temperature. Indeed, the 
latter may well play a role: from our spectral fit, we infer 
Trff = 3139 ± 96 K, while that from its M2 spectral typ e, one 
would estimate T eS = 3514K dLuhman & Riekd [l998). With 
the latter temperature, the object would match the trend much 
better. Because of the large accretion lines in the spectra, the 
models are a poor representation of the data throughout the 
regions selected for our fits. A higher temperature model is 
just a slightly worse fit to the observations. 

5.2. Lithium Line Analysis 

Using the average best-fit effective temperatures and sur- 
face gravities listed in Tables [2]to [71 synthetic spectra, varying 
in lithium abundance from logYVu = 0.0 to 4.0 in steps of 0.5, 
are, for each spectrum, fitted to small spectral regions around 
the 6 104 A and 6 708 A lithium absorption lines. Th e fitting 
method is identical to the method outlined in 34.31 How- 
ever, we constrain the minimization of the least-squares fit 



in a slightly different manner. 

Synthetic spectra of varying lithium abundance are fitted to 
20 A wide spectral regions - specifically 6095-6115 A and 
6695-6715 A. The 6 104 A lithium triplet line is a weaker 
transition than the 6 708 A doublet line. It is also blended 
into the strong 6 103 A Ca II absorption line. As a result, the 
line is only detectable for high lithium abundance. Overall, 
we are not very sensitive to this line, detecting only a small 
change in x 2 ° ver the entire lithium abundance range. This 
is not the case with the stronger lithium doublet at 6 708 A, 
which is very sensitive to lithium abundance, and shows sharp 
transitions from good to bad in its least-squares fits. For high 
lithium abundances, however, the line saturates and without 
the 6 104 A line no good abundance estimates are possible. 

In order to be able to treat all data uniformly, irrespective 
of lithium abundance, we proceeded as follows. First, we use 
the 6 708 A region, with its higher sensitivity, to determine the 
approximate abundance, and select the points with log Vu cor- 
responding to the 4 lowest \ 2 values from the fit to this region. 
For both spectral regions, we then fit a 2nd order polynomial 
to the x 2 f°i" these 4 selected logYVu- Next, we normalize 
both fitted polynomials by dividing by the minimum \ 2 value 
for each spectral region. The two normalized curves are then 
added together to give an average curve, and the minimum 
of this curve is what we take to be the best-fit logYVu- Thus, 
in our procedure for determining the lithium abundances, we 
give equal weight to both regions, unlike the procedure out- 
lined in 34.31 where we desired to keep the weight on the best 
fitted regions by summing the raw x 2 values. 

By way of example, we show, in Fig.|6l the spectral fits to 
the two lithium features for three stars with distinctly different 
lithium abundances. The resulting normalized x 2 curves are 
shown in the bottom panel for each star. Averages of the best- 
fit logYVu for all stars can be found in Tables [2] to [JJ with un- 
certainties representing the error in the mean between multi- 
epoch observations. 

As discussed in 34.51 uncertainties in the model atmo- 
spheres lead to larger errors than the internal errors quoted 
in Tables [2] to [7] Thus, we also handle external errors for the 
lithium abundances with the same approach. We investigate 
changes in the fitted abundances by re-fitting the lithium lines 
to models with r e ff perturbed by ±100 from the initial best-fit 
r en - and also to models with logg also perturbed by ±0.5 dex. 
We find that this parameters leads to changes in logYVu of 
0.15 on average. We show this error on all relevant figures to 
indicate our estimate of the external uncertainties. 

5.3. Ages from the lithium abundance versus temperature 
isochrones 

As a PMS group ages, the lithium abundances of group 
members deplete as a fun ction of luminosity and time (see 
As with EWs in 35.11 we can order the groups in age us- 
ing the relative depletion of members of different groups and 
by comparison with models, we can also determine absolute 
ages. Of course, the absolute ages will only be as good as the 
models. To get an idea of th e as sociated uncertai nty, we try 
both the models of|BCAH9l and|SiessetaLl ((2000J). 

In Fig. [JJ we show the distribution of measured lithium 
abundance as a function of temperature for each of the groups 
in this study, as well as for the field stars (which should be 
fully depleted). As with the EWs in 35.11 it is easy to or- 
der the groups chronologically based on the lithium depletion 
distribution: from oldest to youngest, we again find the order 
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FIG. 6.— Examples of fitting the two lithium features to three stars (from left to right: TWA 02A (M2e), CD-60 416 in TUCHOR (K3/4) and GJ 3305 in BPMG 
(M0. 5)). For each star, the top panel shows the fit to the strong 6 708 A lithium doublet and the middle panel the fit to the much weaker 6 104 A lithium line. 
The bottom panel shows the normalized \ values for the abundances at which we h ad m odels, and the parabolic interpolations used to determine the best-fit 
abundances (which is taken to be at the minimum of the sum of these two curves; see 35 .21 . 
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spectral features as a function of temperature. We identify the ultrafast rota- 
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ABD, TUCHOR, BPMG, and then 77 Cha and TWA both at 
about the same age. It is also clear that the field stars are older 
than all of the PMS groups in this study. 

In Fig. [8] the dependence of logA^i on stellar r e ff is shown 
for ea ch group individually, and a range of PMS is ochrones 
from IBCAH98I (solid lines) and ISiess et all (120001) (dashed 
lines). Both models use a convection mixing length of oimlt = 
1.9 and are scaled to the initial lithium abundance. We chose 
an initial lithium abundance of logyVu = 3.7 t o ma tch the 
abundance we measured from our method (see §15.21 for the 
majority of undepleted stars in the entire sample. The inferred 
ages do depend on this choice, as well as the choice of model. 

From Fig. [8] we can infer ages by eye for each group. We 
tried but decided not to use quantitative analysis, because 
it is apparent that the models do not reproduce the data to 



enough accuracy, especially in the temperature range where 
the lithium abundances are mo st depend ent on age. In gen- 
eral, we find that the models of BCAH98 fit the nature of the 
depletion slightly better for all of the groups, but that the ages 
derived ar e similar, apart from a small systematic offset (with 
IBCAH98I giving slightly higher ages). 

We find that a lithium depletion age of 12±6Myr for 
77 Cha usin g the m odels of IBCAH981 (12±8Myr using the 
ISiess et ail (120001) models; hereafter given in parenthe- 
ses), and and an age of 12±8Myr (12±8M yr) for TWA 
consi s tent with dynamical expansion ages (Jilin ski et al.1 
120051; Ide la Reza et al.l [2006) and other age estimates also 
based onlBCAH98l PMS models dLuhman & Stee"ghl l2004t 
Zuckerman "& Songl 120041 iBarrado y Navascuesl 120061) . For 
BPMG, we find an age of 21±9Myr (13±5Myr), which is 
in ag reement with the estimate of 9-17 Myr from othe r meth- 
ods dacEmanXSongl US |Fei£eiionitaD|2Q06). It is 
slightl y higher than its dy namical expansion age of 11.5 Myr 
from lOrtega et al.l d2002l) . but agrees with a different age 
based on lithium da ting of 10-20 Myr recently found by 
iMamajeketail (120071) . For TUCHOR, we find an age of 
27±llMyr (22±10Myr), which is consistent with all pre- 
vious age estimates for this group (Zuckerman & Webbl l2000T: 
IStelzer & Neuhausedl2000l) . 

For ABD, we find that it is clearly older than TUCHOR 
and clearly younger than the field stars, however the age esti- 
mate from PMS models is poorly constrained. Although the 
field stars show more depletion than ABD, this is not pre- 
dicted by the PMS isochrones; both models show no deple- 
tion after ~ 45 Myr for stars with 4000 < T en < 6000K. Un- 
til the PMS evolutionary models are improved, the best way 
to find an upper limit to the age of the ABD group would 
be to use the cool end of the LDB. This would require stars 
with spectral types later than M3 (T e ff <3 300 K), but, unfor- 
tunately, no such members are known in ABD. Within our 
present large uncertainties, our age estimate is consistent with 
both a y ounger estimate of 50 M yr based on Ha emission 
strength (Zuckerman et al. 2004), as well as an older one of 
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100-140 Myr from H R isochrones dLuhman et al.ll2005l) and 
dynamical expansion (Orteg aet al.ll2007l) l 

We close with a number of notes. First, while the poor fit 
of the data to the models leads to rather large uncertainties on 
the ages, these should be considered overall shifts: the age 
ordering of the groups is secure. Second, an additional un- 
certainty in the derived ages is our choice of initial lithium 
abundance, of logNu = 3.7 based on our observations. De- 
cr easing the initial lithium abundance to logA^i = 3.3 (as used 
in lJeffries & Oliveiral(l2005l) ). predicts younger group ages by 
about 5 Myrs. On the other hand, using a higher initial lithium 
abundance of logNu = 4.0, yields ages larger by 5-10 Myrs. 
Third, we have ignored non-LTE effects in the lithium lines 
(®. While our scatter is larger than the predicted effects, 
the systematic changes with temperature and abundance will 
lead to additional systematic age differences. It also may be 
the underlyi ng reason for our nee d for a relatively high initial 
abundance: ICarlsson et al.l (1 19941) found that around 6000 K, 
the correction for the 6708 line is about -0.3 dex, which would 
imply initial abundances more in line with expectations (at 
these temperatures, the 6 104 A line is very weak and con- 
tributes little to our fits. 

5.4. The Effect of Rotation on Lithium Depletion 

We examine the effect of stellar rotation on lithium 
depletion using projected rotatio nal velocities found pre- 
viously from our observations ( Javaward hana et ail 120061: 
IScholz et al.ll2007h . In Fig. [8] we identify ultrafast rotators as 
those stars that have vsin/ > 70kms _1 . One member (TWA 
6) from TWA, two members (PZ Tel, CD 64dl208) of BPMG 
and three members (CD 53 544, HIP 108422, HIP 2729) of 
TUCHOR are identified as ultrafast rotators. In three of 
the cooler stars (T eS ~3800K), CD64dl208, CD 53 544 and 
HIP 2729 , the lithium EWs are noticeably higher, and the 
derived abundances larger, than the trend in lithium deple- 
tion for the entire group. The correlation between fast ro- 
tation and slower lithium depletion has also be en seen pre- 
vious ly in a sample of weak-line T Tauri stars (Martin et al. 
19941) and in the 115 Myr Pleia des cluster fSoderblo m et al. 
1991 iGarcfa Lopez et alJ[T994b . It may be related to rapidly 
rotating stars bei ng relatively cooler a s the rapid rotation in- 
hibits convection (IChabrier et al.ll2007l) . This alternative view 
is supported by the location of the rotators in Fig. [4] It is evi- 
dent that the T e s derived for these rotators is 1 00-300 K cooler 
than the temperature derived from their spectral types. It may 
be that the presence of a colder equatorial region and a hotter 
polar one affects the model fits differently than the spectral 
typing. 

However, the trend that lithium depletion is slowed down 
by rotation is not seen in all of our ultrafast rotators. The 
relatively slower rotators PZTel (vsin/ = 77.5 km s -1 ) and 
TWA 6 (vsin / = 79.5 km s" 1 ) have lithium abundances com- 
parable to other members in their groups. For our faster ro- 
tator, HIP 108422 (vsin/ =139. 8kms _1 ), the spectrum is so 
strongly broadened that our fitting method does a poor job 
and we see no change in the quality of fit for varying lithium 
abundances. For this, reason, we do not quote an abundance 
for this object, but do note, however, that the equivalent width 
measured is consistent with other group members. 

5.5. Notes on Individual Systems 

BD 17° 6128 - This binary system from BPMG consists of a 
K7 primary with T e ff =4 140 K and logA^Li = 2.78 and a lithium 



depleted secondary of T e ff =3 350 K and \ogNw = 0.54. This 
system is unique to our sample as it is the only case of the cool 
end of the LDB in effect within a binary, and provides a pre- 
cise, if model d ependent, age of the system. Using the models 
from BCAH9"8l an age of 15-50 Myr is predicted, consistent 
with the age we inferred for BPMG as a group. 

GSC 08056-0482 - The measured logjV Li for TUCHOR 
member GSC 08056-0482 is much higher than expected for 
a ~ 30 Myr old, M3 dwarf. Indeed, another M3 star in 
the younger BPMG, GSC 08491-1 194, has an abundance 
over two orders of magnitude smaller (consistent with mod- 
els at ~ 20 Myr). With a modest vsin/ of 34.2 km s" 1 
(I Javawardhana et afll2006l) . rotation does not explain the high 
lithium abundance. As pointed out already b y previous au- 
thors (see Table 3 of Zuckerman & Song 2004), the lithium 
abundance suggests that GSC 08056-0482 is likely not a 
member of TUCHOR, but a star slightly older than TWA and 
r\ Cha, but definitely younger than BPMG. 

6. SUMMARY AND OUTLOOK 

We have measured effective temperatures, surface gravi- 
ties, lithium equivalent widths and lithium abundances for 
121 low-mass PMS stars from five nearby, PMS groups rang- 
ing in age from 8-125 Myr by performing least-squares fits 
of high resolution spectra to sy nthetic spectra created from 
PHOENIX model atmospheres dHauschildt & Baronl [19991) . 
To investigate the reliability of our measurements we com- 
pare the derived r e ff and logg with isochrones from PMS evo- 
lutionary models dBCAH98l) as well as temperatures derived 
from spectral types. 

Isochrones from PMS models for logA^i as a function of 
T e ff are visually compared to the observed distribution. We 
fin d agreeme nt b etween age s deriv ed from PMS isochrones 
of lBCAH98l and ISiess etal] 12000) to ages calculated from 
other methods such as dynamical expansion ages. We find 
that i] Cha and TWA have ages of 12±6Myr and 12±8Myr, 
respectively. BPMG has an age of 21 ±9 Myr, and TUCHOR 
has an age of 27 ±1 1 Myr. We can only constrain a tight lower 
limit for ABD, with an age greater than 45 Myr, since, accord- 
ing to the PMS models, there is no more lithium depletion af- 
ter -45 Myr for stars with 4000 < T eB < 6000K. However, 
the halting of lithium depletion at this age and temperature 
is inconsistent with observations of radial velocity standards 
which demonstrate more depletion than the ABD group and 
the model predictions. Finally, we find that some of the ul- 
trafast rotators in our sample have significantly less lithium 
depletion than other stars in the same group at the same tem- 
perature. 

The consistent determination of T e ff and logg between mul- 
tiple epochs (oj ~ 10 K, u\ ogi — 0.05 dex) means that, in prin- 
ciple, we should be able to constrain those parameters with 
this precision. As revealed by Fig. [3] however, there is an 
apparent systematic offset between the logg inferred from 
model spectra and logg expected from models of stellar evo- 
lution ( 34. 5b . To account for these offsets, we introduce rough 
conservative external errors by examining how the measured 
parameters depend on each other. We find that the systematic 
errors in logg of 0.5 dex lead to systematic errors of 100 K in 
our ability to constrain the T e ff. These external errors, sim- 
ilarly lead to offsets in the measured lithium abundances of 
0.15 dex. 

The small internal errors that we have measured imply that 
currently our accuracy is limited by the models. With further 
improvements in the atmospheric models, there is a potential 



10 





3000 4000 5000 6000 3000 4000 5000 6000 

T ef f (K) T eff (K) 

FIG. 8. — Lithium abundances derived from fits to the spectra as a function of temperature for each group. Overdrawn are predictions from the evolutionary 
models of BCAH98 (solid line) and Siess et al. 12000) (dashed line) as indicated by the ages in the bottom-right corner of each plots. Above the ages, external 
errors are shown. In addition, the ultrafast rotators are identified by black circle outlines, and accretors are identified by black square outlines. 



of comparing r e g and logg directly to evolutionary models, 
thereby finding an age constraint independent of other esti- 
mators, such as color-magnitude diagram or lithium depletion 
boundary fitting. Our data set would be well-suited for use 
with such future improved models. Another use of our data 
set would be to derive both overall metallicity and abundances 
for individual elements. While we do not believe this would 
affect our derived temperatures, etc., to a significant degree, 
it may be interesting to see how uniform the abundances are 
within (and between) groups, and whether there is any depen- 
dence on binarity, etc. 
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TABLE 2 
Results for stars in TW Hydrae 



Object ID 


Sp.T. 


#Obs. 


EW 6708 (A) 


Tes (K) 






logALi 


TWA 1 


K7e 


5 


U.4t>/ ± U.Uzl 


3847 ± 


16 


A AT _L A A"7 
4.0Z ZT U.U / 


3. 


A1 _1_ A A£ 


TWA 2A 


M2e b 


5 


0.567 ± 0.021 


3712 ± 


17 


4.27 ± 0.05 


3 


.36 ± 0.03 


TWA3A 


M3e b 


5 


0.563 ± 0.022 


3268 ± 


12 


3.98 ± 0.05 


3. 


.34 ± 0.05 


TWA 3B 


M3.5 b 


5 


0.529 ± 0.025 


3340 ± 


13 


4.07 ± 0.06 


3. 


.38 ±0.11 


I WA 4AN 




4 


0.359 ± 0.020 


4045 ± 


24 


A C A 1 f\ rv~i 

4.54 ± 0.07 


2 


.40 ± 0.03 


mfA A \ Q 

1 WA 4Ao 




4 


U.44Z ± U.UzU 


4360 ± 


14 


4.ZO ± U.UO 


j. 


1 c i n at 


TWA 5 A 


Ml .5 


5 


U.OU4 ± U.Ozl 


3432 ± 


12 


4.2o ± O.Uj 


3. 


AA _1_ A A/1 

,uy zt 0.04 


TWA 6 


K7 a 


5 


n ^O'i 4-0(171 


3751 ± 


12 


a ah A- OS 


2 


.78 ± 0.04 


TWA 7 


Ml" 


5 


0.546 ± 0.021 


3541 ± 


10 


4.18 ± 0.05 


3 


.42 ± 0.05 


TWA8A 


M2 a 


5 


0.546 ±0.021 


jjj / ztz 


i c 
i j 


4.24 ± 0.05 


3 


.47 ± 0.05 


TWA8B 




5 


0.575 ± 0.026 


3QI3 _j_ 


15 


3.83 ± 0.05 


3 


.54 ±0.11 


TWA9A 


K5 a 


5 


0.535 ± 0.020 


4052 ± 


12 


4.15 ±0.05 


3 


.25 ± 0.03 


TWA9B 


Ml a 


5 


0.517 ±0.021 


3450 ± 


12 


4.19 ±0.05 


3 


.43 ± 0.05 


TWA 10 


M2.5 a 


5 


0.481 ±0.021 


3512 ± 


12 


4.24 ± 0.05 


3 


.34 ± 0.03 


TWA 11B 


M2.5 b 


5 


0.503 ± 0.022 


3653 ± 


14 


4.26 ± 0.05 


3 


.22 ± 0.03 


TWA 12 


M2 b 


5 


0.521 ± 0.020 


1{\A1 4- 

JOT / Hi 


I t. 


4.33 ± 0.05 


3. 


.17 ±0.03 


TWA 13A 


Mle b 


5 


0.533 ± 0.021 


3845 ± 


Lo 


4.46 ± 0.05 


3 


.23 ± 0.03 


TWA 13B 


M2e b 


5 


0.577 ± 0.021 


3817 ± 


16 


4.47 ± 0.05 


3. 


.36 ± 0.03 


TWA 14 


M0 b 


5 


0.589 ± 0.021 


3701 ± 


13 


4.45 ± 0.05 


3 


.63 ± 0.06 


TWA 1 S A 
1 vva ija 


Ml 5 b 




U.t7 i r ZL U.UZJ 


3482 ± 


12 


A -|_ A OS 


3 


07 -t- A OS 


1 WA 1 JO 


M2 b 


5 


A AQ A _L n ATA 


3483 ± 


12 


A 1A _L A A^ 

4. JU ± U.UO 


3 


1 £ _1_ A AC 

. lo zt 0.05 


TWA 16 


M1.5 b 


5 


0.354 ± 0.024 


3610 ± 


20 


4.27 ± 0.05 


2 


.55 ± 0.04 


TWA 17 


K5 b 


5 


0.499 ± 0.021 


3884 ± 


25 


4.02 ± 0.06 


3. 


.31 ±0.10 


TWA 18 


M0.5 b 


5 


0.464 ± 0.021 


3780 ± 


13 


4.50 ± 0.06 


3. 


.14 ±0.04 


TWA 19A 


G5 a 


5 


0.191 ±0.020 


5986 ± 


14 


3.83 ± 0.07 


3 


.92 ± 0.03 


TWA 19B 


K7 b 


5 


0.452 ± 0.022 


3896 ± 


15 


4.24 ± 0.05 


2 


.89 ± 0.07 


TWA 21 


K3/4 a 


5 


0.369 ± 0.020 


4689 ± 


14 


4.44 ± 0.05 


3. 


.26 ± 0.03 


TWA 22 


M5 a 


5 


0.616 ±0.021 


2990 ± 


13 


4.20 ± 0.05 


3 


.57 ± 0.06 


TWA 23 


Ml a 


5 


0.525 ± 0.020 


3466 ± 


12 


4.12 ±0.05 


3. 


.51 ±0.04 


TWA 24N 




5 


0.438 ± 0.032 


3669 ± 


15 


4.28 ± 0.06 


2 


.83 ± 0.08 


TWA 24S 


K3 a 


5 


0.380 ± 0.020 


4920 ± 


16 


4.32 ± 0.05 


3 


.64 ± 0.02 


TWA 25 


M0 a 


1 


0.570 ± 0.020 


3920 ± 


10 


4.45 ± 0.05 


3 


.25 ± 0.02 



NOTE. — all uncertainties are internal, derived from the scatter of fitted values 
from individual spectra from the mean. The external uncertainties are much larger 

( see &31>. 

a Zuckerman & Sond j2004l) b lde la Reza et all d2006h 



TABLE 3 

Results for stars in r\ Chamaeleontis 



Object ID 


Sp.T. 


#Obs. 


EW 6708 (A) 


T ets (K) 


logs 


logMj 


ryCha 1 


K4 a 


5 


0.511 ±0.020 


4107 ± 15 


3.94 ± 0.05 


3.31 ± 0.05 


■q Cha 3 


M3.25 b 


5 


0.542 ± 0.021 


3508 ± 12 


4.18 ±0.05 


3.55 ± 0.04 


?jCha4 


K7 a 


5 


0.587 ± 0.021 


3822 ± 1 1 


4.43 ± 0.05 


3.32 ± 0.03 


r\ Cha 5 


M4 b 


5 


0.606 ±0.021 


3314 ± 11 


4.04 ± 0.05 


3.41 ± 0.09 


r\ Cha 6 


M2 a 


4 


0.489 ± 0.022 


3492 ± 1 1 


4.22 ± 0.06 


3.31 ±0.03 


?7Cha7 


K6 b 


5 


0.415 ± 0.022 


4002 ± 25 


3.81 ± 0.05 


2.77 ± 0.08 


r?Cha9 


M4.5 a 


5 


0.566 ± 0.027 


3127 ± 11 


3.84 ± 0.05 


3.36 ±0.11 


jjCha 10 


K7 a 


5 


0.534 ±0.021 


3992 ± 11 


4.67 ± 0.06 


3.08 ± 0.03 


T) Cha 11 


K4 a 


5 


0.483 ±0.021 


4182 ± 19 


3.92 ± 0.05 


3.28 ± 0.03 


r?Cha 12 


M2 a 


5 


0.609 ± 0.023 


3440 ± 11 


4.31 ± 0.06 


3.20 ± 0.08 


r?Cha 13 


M2 a 


5 


0.426 ± 0.023 


3139 ± 11 


3.68 ± 0.05 


2.89 ± 0.07 



NOTE. — all uncertainties are internal, derived from the scatter of fitted values 
from individual spectra from the mean. The external uncertainties are much larger 

(seej43L 

a Zuckerma n & Sond d2004h b lLuhrrian & Steeghs d2004h 



TABLE 4 

Results for stars in the (3 Pictoris moving group (BPMG) 



Ohiprt JT) 
(JUJCL.I 1 iJ 


Sp T 


# Obs 


11 "6708 K A ) 


Tiff (K) 






AO Men 


K6/7 


5 


0.409 ± 0.020 


4317 ± 


12 


4.13 ±0.05 


3.01 ± 0.02 


AUMic 


Mle 


5 


0.067 ± 0.020 


3841 ± 


1 1 


4.66 ± 0.05 


0.80 ± 0.03 


BD-17 6128A 


K7e/M0 


5 


0.427 ± 0.020 


4140 ± 


1 3 


4.37 ± 0.05 


2.78 ± 0.02 


BD-17 6128B 




4 


0.051 ± 0.035 


3350 ± 


i n 
19 


4.28 ± 0.05 


0.54 ±0.10 


CD -64 1208 


K7 


5 


0.464 ± 0.021 


3740 ± 


19 


3.96 ± 0.05 


1.27 ± 0.04 


GJ 3305 


M0.5 


5 


0.077 ± 0.020 


3840 ± 


1 1 


4.72 ± 0.05 


1.09 ±0.03 


GJ 799N 


M4.5e 


5 


0.020 ± 0.020 


3237 ± 


10 


4.47 ± 0.05 


0.25 ± 0.02 


GJ 799S 


M4.5e 


5 


0.020 ± 0.021 


3244 ± 


1 1 


4.46 ± 0.05 


0.25 ± 0.02 


HD 164249B 




4 


0.075 ± 0.029 


3505 ± 


1 8 


4.32 ±0.10 


0.25 ±0.10 


HD 181327 


F5.5 


5 


0.122 ± 0.020 


6392 ± 


1 1 


4.17 ±0.08 


3.70 ± 0.02 


HD 35850 


F7 


5 


0.146 ± 0.020 


6275 ± 


18 


4.15 ±0.08 


3.68 ± 0.02 


HIP 10679 


G2V 


1 


0.172 ±0.020 


6080 ± 


10 


4.25 ± 0.05 


3.70 ± 0.02 


HIP 10680 


F5V 


1 


0.132 ±0.020 


6250 ± 


10 


4.00 ± 0.05 


3.65 ± 0.02 


HIP 112312A 


M4e 


3 


0.020 ± 0.020 


3340 ± 


12 


4.37 ± 0.07 


0.25 ± 0.02 


HIP 112312B 


M4.5 


3 


0.315 ±0.022 


3170 ± 


14 


4.17 ±0.07 


2.17 ±0.19 


HIP 11437A 


K8 


1 


0.248 ± 0.020 


4310 ± 


10 


4.60 ± 0.05 


2.30 ± 0.02 


HIP 11437B 


M0 


1 


0.130 ±0.020 


3620 ± 


10 


4.55 ± 0.05 


1.30 ±0.02 


HIP 12545 


M0 


2 


0.433 ± 0.023 


4205 ± 


11 


4.22 ± 0.06 


2.92 ± 0.03 


HIP 23309 


M5 


5 


0.336 ± 0.022 


4041 ± 


10 


4.72 ± 0.05 


2.46 ± 0.02 


HIP23418N 


M3V 


2 


0.020 ± 0.021 


3335 ± 


11 


4.43 ± 0.06 


0.25 ± 0.02 


HIP23418S 




2 


0.020 ± 0.022 


3310 ± 


14 


4.35 ± 0.07 


0.25 ± 0.02 


PZTel 


KOVp 


5 


0.271 ±0.020 


5537 ± 


15 


3.98 ±0.11 


3.65 ± 0.02 


V343 Nor 


K0V 


5 


0.286 ± 0.020 


5520 ± 


13 


4.22 ± 0.07 


3.75 ± 0.02 



NOTE. — all uncertainties are internal, derived from the scatter of fitted values from 
individual spectra fro m the mean. The exte r nal un certainties are much larger (see 34.5^ 
All spectral types from Zuckerma n & Sonpl {2004) 
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TABLE 5 

Results for stars in Tucanae-Horologium (TUCHOR) 



Object ID 


Sp.T. 


#Obs. 


EW 6708 (A) 


T ef f (K) 




logMj 


/T"\ CI CA A 

L.U -jj 544 


KoVe 


3 


U.2 /;> ± U.U21 


3796 ± 


13 


4.33 ± 0.05 


l.ou it 0.04 


CJJ -60 416 


VI IA 

ro/4 


c 

J 


n at"7 _l_ A aa 1 
U.ZJ / ± U.U21 


43 1 ± 


33 


4.34 it O.Oo 


Z.33 it 0.03 


/TJTTi A/1 1 On 

CrD -64 1 20 


Kl Ve 


5 


a ao/1 J_ n mi 
U.254 ± U.U2T 


5212 ± 


48 


4.ol ± 0.03 


Q HI _L A f\1 

j. 15 ± 0.03 


r^cr^ Qnc£ ri^ioi 
LriSL. o056-04o2 


M3Ve 


3 


U.Jjo ± U.U21 


3541 ± 


13 


a K.A _L n r\z 
4.34 it 0.03 


Z.38 it 0.06 


LrSL. 8491-1 194 


M3Ve 


5 


A AT 1 1 A AT 1 

U.U3J ± 0.021 


3578 ± 


12 


4.64 ± O.Oo 


0.25 ± 0.02 


/■ ' i 1 a ^ o ac\~i annc 

LrSL. 849 /-0995 


K6Ve 


4 


A 1 AT 1 A ATA 

U.127 ± 0.022 


4238 ± 


1 8 


4.64 ± 0.0/ 


1 CI 1 A AC 

1.52 ± 0.05 


Tjri 1 1 1 OQ 

rilJ 131o3 


G5V 


5 


A AA£ J_ A AAA 
U.2Uj it U.U2U 


5854 ± 


22 


4.Z0 it 0.03 


1 HI _L A AO 

3. id ± 0.02 


rilJ 13246 


F8V 


c 

J 


(1 (]] _1_ A AAA 

U. Ijj ± u.uzu 


6292 ± 


13 


A 1 1 _l_ A AA 

4. 1 / it 0.06 


1 AO _L A (Yl 

3.65 ± 0.03 


rilJ ojjo 


AW 

Ltd V 


c 

J 


A 1 OA _1_ A AAA 

u. iyz ± u.uzu 


5825 ± 


19 


4.Z / it 0.03 


i oa J_ n no 
3.90 it 0.02 


T TT~\ AAC A 

tiu 9054 


K1V 


5 


A 1 TO 1 A AAA 

U.l/o ± U.U2U 


5165 ± 


35 


1 /r 1 A AC 

4.65 it 0.05 


3.01 ± 0.02 


T TTT1 1 A C "f O O 

HIP 105388 


G5V 


5 


A AA A 1 A AAA 

0.224 ± 0.020 


5748 ± 


1 1 


a 1 o I n no 

4. 18 it 0.08 


3.71 ± 0.02 


TJTD 1 AC A A/1 

Hlr 105404 


K0V 


3 


A 1 "7 1 _L_ A AAA 
U. 1 / 1 ± U.U22 


5543 ± 


54 


A TO _L A 1 A 

4.Z5 it 0.10 


1 CA _L A AC 

3.56 ± 0.05 


tjtd 1 nil A c 

Hlr 10/345 


M 1 


5 


A ACC J_ A AA 1 
U.UjJ ± U.U21 


3975 ± 


13 


A OO _L A ri5 

4.oz it 0.03 


A O C _L A AA 

0.25 it 0.09 


tjtd iacioo 
rilr 10o422 


G8V 




A A/^f _L- A AAA 
U.2ol ± U.U2U 


5541 ± 


27 


1 AO _L A AT 

3.99 ± 0.0/ 


N/A a 


T TTT) 1 1 1 1 

HIP 1 1 1 3 


G6V 


3 


A AAA 1 A AAA 

0.202 ± 0.020 


5757 ± 


13 


4.35 ± 0.0/ 


Ton 1 f\ l\A 

3.80 ± 0.04 


T_TT~D 1 1 AO/1 QM 

Hlr 1 16/4oJN 




3 


A A 1 Q _1_ A AAA 

U.21o ± U.UzU 


4620 ± 


10 


A AH _L A f\Z 

4.4/ it 0.05 


1 AO _L A fV2 

2.62 it 0.03 


UTTJ 1 1 A "7 /IOC 

rilr 1 1 6/49.3 


Lr51 V 


i 
j 


A A 1 A -1- A AA 1 

u.ziz ± u.uzi 


5813 ± 


44 


A in _U A A A 
4.30 it 0.06 


A AA J_ A AO 

4.00 it 0.02 


IJTTJ 1/10 1 

rilr 14ol 


CO 

ro 


3 


A 1 TO _L- A AAA 

U. 12o ± U.UZU 


6323 ± 


13 


A T7 _L A AA 
4.2/ it 0.06 


1 AO _L A A1 

3.6o ± 0.03 


Hlr 16853 


G2V 


5 


A 1 /I A 1 A ATA 

0. 149 ± 0.020 


6217 ± 


14 


4.29 ± 0.07 


A An 1 r\ AO 

4.00 ± 0.02 


tjtd mm 
rilr lylU 


M 1 


3 


A 1 O 1 J_ A ATA 

U. lol ± U.U2U 


3890 ± 


20 


J 71 1 A A^ 

4. / J ± U.06 


1 AT l A A^ 


rilr 1993 


Ml 


3 


A (110 A AAT 

U.Uio ± U.U22 


4017 ± 


16 


/I OT ^ A AO 

4.0T ± U.Uo 


A A^ -L- A AA 
U.2j it U.U2 


rilr 21632 


G3V 


5 


A ICQ J_ A ATA 

U. loo it U.UzU 




1 8 


A 1 T _L A A^ 

4. 12 it U.Uo 


1 OT _L A AA 
J.02 ± U.U2 


rilr 22295 


F7V 


5 


A 1 TA _L A ATA 
U. ItU it U.U2U 


6342 ± 


17 


/I 1 1 _1_ A AO 

4. i j it u.uy 


Q ^ ^ J_ AAT 
J.DJ it U.U2 


rilr 2/29 


K5V 


3 


A 1 1 Q J— A ATA 

U.Jto it U.UzU 


3827 ± 


11 


/l AA _L A A£ 
4.U2 it U.Uo 


A ^O J— A AT 

U.DO it U.Ut 


HIP 30030 


GO 


5 


0.163 ±0.021 


6210 ± 


21 


4.19 ± 0.08 


3.77 ± 0.04 


HIP 30034 


K2V 


5 


0.287 ± 0.020 


5268 ± 


23 


4.69 ± 0.06 


3.57 ± 0.02 


HIP 32235 


G6V 


5 


0.233 ± 0.020 


5774 ± 


13 


4.34 ± 0.07 


3.99 ± 0.02 


HIP 33737 


K3V 


5 


0.279 ± 0.020 


4859 ± 


18 


4.64 ± 0.05 


3.16 ±0.03 


HIP 3556 


M3 


3 


0.055 ± 0.022 


3677 ± 


17 


4.72 ± 0.05 


0.27 ± 0.09 


HIP 490 


GOV 


3 


0.153 ±0.020 


6173 ± 


17 


4.28 ± 0.05 


4.00 ± 0.02 


HIP 9141 


G3/5V 


5 


0.187 ±0.020 


5992 ± 


13 


4.23 ± 0.05 


3.80 ± 0.02 


TUCh 7600-0516 


Kl 


5 


0.249 ± 0.020 


5163 ± 


31 


4.63 ± 0.06 


3.56 ± 0.02 


TYC 5882-1 169 


K3/4 


5 


0.241 ± 0.022 


4939 ± 


34 


4.65 ± 0.09 


3.05 ±0.11 


TYC 7065-0879N 


K4 


3 


0.259 ± 0.020 


5513 ± 


28 


4.58 ± 0.07 


3.63 ± 0.03 


TYC 7065-0879S 


K4 


3 


0.267 ± 0.021 


5453 ± 


31 


4.70 ±0.14 


3.63 ± 0.03 



NOTE. — all uncertainties are internal, derived from the scatter of fitted values from 
individual spectra fro m the mean. The external u ncertainties are much larger (see 94.51 1. 
All spectral types from Zuckerma n & Songl (2004) 

a Lithium abundance could not be measured due to high rotational broadening in 
HIP 108422 (vsini = 139.80 km/s) 
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TABLE 6 

Results for stars in AB Doradus (ABD) 



Ohif rt IT) 


jp. l . 


# Obs 


^"6708 l/V 


Tat (K) 


l°g£ 




AB Dor 


Kl 


3 


0.261 ± 0.021 


5210 ± 


34 


4.63 ± 0.06 


3.50 ± 0.04 


GSC 08894-00426 


M2 


3 


0.037 ± 0.022 


3343 ± 


1 1 


4.83 ± 0.05 


0.25 ± 0.02 


HD 13482A 


Kl 


2 


0.124 ±0.020 


5725 ± 


105 


4.47 ±0.18 


3.17 ±0.03 


HD 13482B 




1 


0.076 ± 0.020 


5090 ± 


10 


4.95 ± 0.05 


1.75 ± 0.02 


HD 17332B 




1 


0.168 ±0.020 


5780 ± 


10 


4.35 ± 0.05 


3.55 ± 0.02 


HD 217343 


G3V 


3 


0.172 ±0.020 


5957 ± 


28 


4.25 ± 0.05 


3.68 ± 0.04 


HD 217379N 




2 


0.020 ± 0.020 


4050 ± 


14 


5.07 ± 0.06 


0.30 ± 0.05 


HD217379S 




2 


0.030 ± 0.021 


4165 ± 


1 1 


5.05 ± 0.07 


0.35 ± 0.05 


HD 218860 


G5 


2 


0.222 ± 0.020 


5670 ± 


10 


4.32 ± 0.06 


3.55 ± 0.02 


HD 224228 


K3V 


2 


0.076 ± 0.020 


4905 ± 


11 


4.82 ± 0.06 


1.35 ± 0.02 


HD 35650 


K7 


4 


0.020 ± 0.020 


4252 ± 


10 


5.06 ± 0.05 


0.25 ± 0.02 


HD 45270 


GIV 


3 


0.142 ±0.020 


6187 ± 


24 


4.50 ±0.10 


3.90 ± 0.04 


HD 65569 


KIV 


3 


0.155 ±0.020 


5223 ± 


31 


4.87 ± 0.05 


2.77 ± 0.03 


HIP 14807 


K6 


1 


0.034 ± 0.020 


4500 ± 


10 


4.60 ± 0.05 


0.55 ± 0.02 


HIP 14809 


G5 


1 


0.150 ±0.020 


6050 ± 


10 


4.25 ± 0.05 


3.60 ± 0.02 


HIP 17695 


M3 


2 


0.074 ± 0.029 


3545 ± 


18 


4.72 ± 0.06 


0.25 ± 0.02 


HIP 26369 


K7 


3 


0.044 ± 0.020 


4010 ± 


15 


4.55 ± 0.06 


0.25 ± 0.02 


HIP 31878 


K7 


3 


0.042 ± 0.024 


4240 ± 


14 


5.12 ±0.05 


0.53 ± 0.03 


HIP 6276 


G8 


3 


0.153 ±0.020 


5643 ± 


13 


4.43 ± 0.05 


3.22 ± 0.04 


HR 2468 


G1.5 


2 


0.138 ±0.020 


6060 ± 


14 


4.60 ±0.11 


3.58 ± 0.03 


UYPic 


KOV 


3 


0.267 ± 0.020 


5600 ± 


15 


5.05 ± 0.06 


3.60 ± 0.02 


V372 Pub 


M3 


5 


0.020 ± 0.021 


3688 ± 


19 


4.88 ± 0.05 


0.25 ± 0.02 



NOTE. — all uncertainties are internal, derived from the scatter of fitted values from 
individual spectra from the mean. The external u ncertainties are much larger (see 34. 5\ . 
All spectral types from Zuckerman & Songl {2004) 



TABLE 7 

Results for the field stars (radial velocity standards) 



Object ID 


Sp.T. 


#Obs. 


7kf (K) 




log/Vu 


GJ729 


M3.5 


1 


3340 ± 10 


5.00 ± 0.05 


0.25 ± 0.02 


GJ 156 


K7 


4 


4110 ± 11 


5.19 ±0.05 


0.25 ± 0.02 


Gl 205 


M1.5 


4 


3875 ± 10 


4.89 ± 0.06 


0.25 ± 0.02 


Gl 349 


K3 


5 


4712 ±24 


4.79 ± 0.07 


0.35 ± 0.02 


G1382 


M1.5 


2 


3745 ± 11 


4.68 ± 0.06 


0.25 ± 0.03 


Gl 876 


M4 


1 


3290 ± 10 


4.60 ± 0.05 


0.25 ± 0.02 


G1880 


M1.5 


2 


3795 ± 11 


4.97 ± 0.06 


0.25 ± 0.02 


HD 103932 


K5 


5 


4360 ± 10 


4.77 ± 0.05 


0.25 ± 0.02 


HD 111631 


K7 


5 


4073 ± 10 


5.19 ±0.05 


0.25 ± 0.02 


HD 153458 


GO 


2 


5935 ± 1 1 


4.12 ±0.06 


2.25 ± 0.02 


HD 172051 


G5V 


3 


6070 ± 12 


4.45 ± 0.06 


0.97 ± 0.03 


HD 120467 


K4 


5 


4231 ± 11 


4.97 ± 0.05 


0.25 ± 0.02 


HD 83443 


K0 


3 


5280 ± 12 


3.90 ± 0.06 


0.45 ± 0.02 


HD 87359 


G5 


3 


5903 ± 35 


4.15 ±0.05 


0.92 ± 0.04 


HD 88218 


F8 


1 


6240 ± 10 


4.10 ±0.05 


2.80 ± 0.02 


HD 92945 


Kl 


1 


5310 ± 10 


4.70 ± 0.05 


2.95 ± 0.02 


HD 96700 


G2 


1 


6160 ± 10 


4.20 ± 0.05 


1.55 ±0.02 


LHS 1763 


K5 


4 


4348 ± 10 


4.78 ± 0.05 


0.49 ± 0.02 


NSV 2863 


M1.5 


4 


3882 ± 11 


4.94 ± 0.06 


0.25 ± 0.02 


NSV6431 


M2 


5 


3724 ± 10 


4.64 ± 0.05 


0.25 ± 0.02 



NOTE. — all uncertainties are internal, derived from the scatter 
of fitted values from individual spectra from the mean. The external 
uncertainties are much larger (see 34.5b . 



